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00 1 ABSTRACT 

(N ' 

Normal and superconducting state spectral properties of cuprates are theoretically described within the ex- 
tended t-J model. The method is based on the equations of motion for projected fermionic operators and the 
mode-coupling approximation for the self-energy matrix. The dynamical spin susceptibility at various doping 
is considered as an input, extracted from experiments. The analysis shows that the onset of superconductivity 
is dominated by the spin-fluctuation contribution. The coupling to spin fluctuations directly involves the next- 
nearest- neighbor hopping t', hence T c shows a pronounced dependence on t' . The latter can offer an explanation 
for the variation of T c among different families of hole-doped cuprates. A formula for maximum T c is given and 
it is shown that optimum doping, where maximum T c is reached, is with increasing — t' progresively increased. 
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; 1. INTRODUCTION 

The mechanism of high-temperature supeconductivity (SC) in cuprates represents one of the central open ques- 
tions in the solid state theory. The role of strong correlations and the antifcrromagnetic (AFM) state of the 
reference insulating undoped compound has been recognized very early. 1 Still, up to date there is no general con- 
sensus whether ingredients as embodied within the prototype single-band models of strongly correlated electrons 
, are sufficient to explain the onset of high T c , or in addition other degrees of freedom, as e.g. phonons, should be 
1 invoked. As the basis of our study, we assume the extended t-J model, 2 allowing for the next-nearest-neighbor 
(NNN) hopping t' term. The latter model, as well as the Hubbard model, 4 both closely related in the strong 
correlation limit U 3> t, have been considered by numerous authors to address the existence of SC due to strong 
correlations alone. Within the parent resonating-valence-bond (RVB) theory 1-3 and slave-boson approaches to 
' the t-J model 5 the SC emerges due to the condensation of singlet pairs, induced by the exchange interaction J. 
An alternative view on strong correlations has been that AFM spin fluctuations, becoming particularly longer- 
ranged and soft at low hole doping, represent the relevant low-energy bosonic excitations mediating the attractive 
interaction between quasiparticles (QP) and induce the d-wave SC pairing. The latter scenario has been mainly 
followed in the planar Hubbard model 4 and in the phcnomcnological spin-fermion model. 6 Recent numerical 
^ , studies of the planar t-J model using the variational quantum Monte Carlo approach, 7 as well as of the Hubbard 
model using cluster dynamical mean-field approximation, 8 seem to confirm the stability of the d-wave SC as the 
ground state at intermediate hole doping. The relevance of t' for T c has been already recognized 9 and recently, 
O ■ there are also some numerical studies of the influence of t' on pairing within prototype models, 10:11 although 
, with conflicting conclusions. 

^ ■ 2. MODEL AND METHOD 

3 • 

The t-J model is nonperturbative by construction, so it is hard to design for it a trustworthy analytical method. 
One approach is to use the method of equations of motion (EQM) to derive an effective coupling between 
fermionic QP and spin fluctuations. 12 The latter method has been employed to evaluate the self energy and 
anomalous properties of the spectral function, 12-14 in particular the appearance of the pseudogap and the 
effective truncation of the Fermi surface (FS) at low hole doping. 14 The analysis has been extended to the study 
of the SC pairing, 13 ' 15 while an analogous approach has been also applied to the Hubbard model. 16 
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2.1. Equations of motion 

We employ the formalism of the EQM and the resulting Eliashbcrg equations within the simplest mode-coupling 
approximation. 13-15 Equations involve the dynamical spin susceptibility which we consider as an phenomcnolog- 
ical input taken from the inelastic- neutron-scattering (INS) and NMR- relaxation experiments in cuprates. The 
analysis of these experiments 17 reveals that in the metallic state the AFM staggered susceptibility is strongly 
enhanced at the crossover from the overdoped (OD) regime to optimum (OP) doping and is increasing further in 
underdoped (UD) cuprates, while at the same time the corresponding spin-fluctuation energy scale is becoming 
very soft. Direct evidence for the latter is the appearance of the resonant magnetic mode 18,19 within the SC 
phase indicating that the AFM paramagnon mode can become even lower than the SC gap. These facts give 
a support to the scenario that spin fluctuations in cuprates represent the lowest bosonic mode relevant for the 
d-wave SC pairing. 

One of the central results of our EQM approach is that the relevant coupling to AFM paramagnons involves 
directly t' , but not t. The evident consequence is the sensitivity of T c on t', consistent with the experimental 
evidence for different families of cuprates. 20 

We consider the extended t-J model 



including both the NN hopping tij — t and the NNN hopping tij = t'. The projection in fermionic operators, 
c is = (1 — rii- s )ci S leads to a nontrivial EQM, which can be in the k basis represented as 

[c k5 , H] = [(1 + c h )^ - J(l - c h )]c ks + (2) 



where Ch is the hole concentration and mk q is the effective spin-fermion coupling mk q = 2J7 q + e k _ q , while 
e k = — 4f7 k — 4i'7 k is the bare band dispersion on a square lattice. We use the symmetrized coupling as derived 
in Rcf. 14 to keep a similarity with the spin-fermion phenomenology 6 

m k q = 2J7q + ^(e k _q + e k ). (3) 

EQM, Eq. (2), are used to derive the approximation for the Green's function (GF) matrix G k s(w) = ((^ksl^ks))" 
for the spinor = (ck, s , c^ k _ s ). 

2.2. Gap equation 

We follow the method, as applied to the normal state (NS) GF by present authors, 12, 14 and generalized to the 
SC pairing in Ref. 13. In general, we can represent the GF matrix in the form 

Gk s (w) _1 = ^[^to - Cks + M r 3 - EksM], ( 4 ) 

where a = 5Zj({cj S , c| } + )/iV = (1 + c/j)/2 is the normalization factor, [i is the chemical potential and the 
frequency matrix, (ks = ^({['I'ks, H], ^t_ s }+), which generates a renormalized band C k = Ctl = C ~ ^VitDi — 
4??2t'7k an d the mean-field (MF) SC gap 

K = C k ? - --^^7k- q (c- q ,- s c q;S ). (5) 
q 

To evaluate E ks (w) we use the lowest-order mode-coupling approximation, analogous to the treatment of the 
SC in the spin-fermion model, 6 introduced in the t-J model for the NS GF 12, 14 and extended to the analysis of 



the SC state. 13 Taking into account EQM, Eq. (2), and by decoupling fermionic and bosonic degrees of freedom, 
one gets 

^11(12)/. n —3 v ^ 2 ,-,11(12)/. \ /• \ f C \ 

S ks ( lUJ n) = 2^ m kq G k-q,s(^m)Xq(^„ - ^m) (6) 

q,m 

where iuu n = iir(2n+ l)//3 and Xq( w ) is t ne dynamical spin susceptibility, whereby we have neglected the charge- 
fluctuation contribution. 

In order to analyze the low-energy behavior in the NS and in the SC state, we use the QP approximation for 
the spectral function matrix 

^ksM ~ ^^(wto - e k T 3 - A ks ri)[5(w - Ek) - 5(u) + E k )], (7) 
Zhik 

where E k = (e k + A ks ) x / 2 , while NS parameters, i.e., the QP weight Z k and the QP energy ek, are determined 
from G ks (uj ~ 0), Eq. (4). The renormalized SC gap is 

A ks = Z k [A k + £ 12 (0)]. (8) 

It follows from Eq. (4) that G k 2 (iu) n ) ~ —aZ^A^s/ (uj 2 + E 2 .). By defining the normalized frequency dependence 
Fq(iu)i) — Xq(* w 0/Xqi an d rewriting the MF gap, Eq. (5), in terms of the spectral function, Eq. (7), we can 
display the gap equation in a more familiar form, 



iks = ^H[ 4J7k -1~ 3r "k : k-qXk-qCq,k- 



q 

(Z°Z°A qs /2i? q )th(/?iV2), (9) 
where C kq = 7 kq (iw„ ~ 0)/J k plays the role of the cutoff function with 

/kq(iw„) = i F ^ iuJ n ~ iUm) — 2 , (10) 

m m ks 

and J k = th(/3_E k /2)/(2_E k ). Eq. (9) represents the BCS-like expression which we use furtheron to evaluate T c , 
as well as to discuss the SC gap A q (T = 0). To proceed we need the input of two kinds: a) the dynamical spin 
susceptibility Xq( w )> an d b) the NS QP properties Z k ,ek. 

2.3. Parameters 

The INS experiments show that within the NS the low-w spin dynamics at q <~ Q is generally overdamped in 
the whole doping (but paramagnetic) regime. 19 Hence we assume Xq( w ) °f t ne form 

*» " ^TJ' *■«■*> " Rfe- <"> 

Following the recent memory-function analysis 21 £? q = x q Fq should be quite independent of q = q — Q. We 
choose the variation as r q ~ Fq(1 + wq 2 / k 2 ) 2 consistent with the INS observation of faster than Lorentzian 
fall-off of X q ( w ) vs - Q- 19 w ~ 0-42 in order that k represents the usual inverse AFM correlation length. 

Consequently, we end up with parameters Xq^q,k, which are dependent on Ch, but in general as well 
vary with T. Although one can attempt to calculate them using the analogous framework, 21 we use here the 
experimental input for cuprates. We refer to results of the recent analysis, 17 where NMR T2G relaxation and 
INS data were used to extract k, XqCO an d ^q(T) for various cuprates, ranging from the UD to the OD regime. 
For comparison with the t-J model, we use usual parameters t = 400 meV, J = 0.3t. At least for UD cuprates, 
quite consistent estimates for Xq,Tq can b e obtained also directly from the INS spectra. 19 For UD, OP and 
OD regime, i.e., Ch = 0.12,0.17,0.22, respectively, we use furtheron the following values: Xq^ = 15.0,4.0,1.0, 
T° Q /t = 0.03,0.1,0.18 (appropriate at low T), and k = 0.5,1.0,1.2. It is evident, that in the UD regime the 
energy scale Tq becomes very small (and consequently Xq 1/Tq targe, in spite of modest k 17 ), supported by 
a pronounced resonance mode. 19 We take into account also the T dependence, i.e., Tq(T) ~ Tq + T, 17 being 
significant only in the UD regime. 
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Figure 1. QP weight evaluated for t'/t = —0.3 for parameters corresponding to c/, = 0.12 and = 0.22, respectively. 
White line represents the location of the FS. 



3.1. Normal state 

For the NS Ak{uS) and corresponding Zk,£k we s °l ve Eq. (6) for E^ 1 = £k as in Ref. 14, with the input for 
Xq(w) as described above. Since our present aim is on the mechanism of the SC, we do not perform the full 
self-consistent calculation of Sk(w), but rather simplify it as done in the previous study. 14 Large incoherent 
Ek(w <C 0) leads to an overall decrease of the QP weight Z < 1 and the QP dispersion with renormalized 
771 , 772 < 1, which we assume here following Ref. 14, as 771 = r\i = 0.5, Z — 0.7. Soft AFM fluctuations with 
q ~ Q lead through Eq. (6) to an additional reduction of Z^, which is k-dependent. A pseudogap appears along 
the AFM zone boundary and the FS is effectively truncated in the UD regime with Zk F <C 1 near the saddle 
points (tt, Q) (in the antinodal part of the FS). 14 We fix ^ with the FS volume corresponding to band filling 



The coupling to low-energy AFM fluctuations, Eq. (6), leads to an additional QP renormalization. For fixed 
t'/t = —0.3, we present in Fig. 1 results for the variation of the Z^ in the Brillouin zone for two sets of parameters, 
representing the UD and the OD regime, respectively. The location of the renormalized FS is also presented in 
Fig. 1. While the coupling to AFM fluctuations partly changes the shape of the FS, more pronounced effect is on 
the QP weight. It is evident from Fig. 1 that Z^ is reduced along the AFM zone boundary away from the nodal 
points. Particularly strong renormalization < 1 happens in the UD case, leading to an effective truncation 
of the FS away from nodal points. 14 

3.2. Superconducting state 

First we comment on the general properties of the gap equation, Eq. (9). Close to half-filling and for % q peaked 
at q ~ Q both terms favor the d x 2_ y 2 SC. The MF-part A£, Eq. (5), involves only J which induces a nonretarded 
local attraction, playing the major role in the RVB theories. 1,2 In contrast, the spin-fluctuation part represents 
a retarded interaction due to the cutoff function Ck q determined by Tk-q. The largest contribution to the SC 
pairing naturally arises from the antinodal part of the FS. Meanwhile, in the same region of the FS also Zk is 
smallest, reducing the pairing strength in particular in the UD regime. Our analysis is also based on the lowest 
order mode-coupling treatment of the SC pairing as well as of the QP properties near the FS. Taking this into 
account, one can question the relative role of the hopping parameters t, t' and the exchange J in the coupling, 
Eq. (3). While our derivation within the t-J model is straightforward, an analogous analysis within the Hubbard 



3. NUMERICAL RESULTS 
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Figure 2. T c /t vs. doping c h for t'/t = —0.3, calculated for various versions of Eq. (9): a) full result (full line), b) with 
neglected MF term (dashed line), and c) in addition to b) modified m kq without the J term (dotted line). 

model using the projections to the lower and the upper Hubbard band, respectively, would not yield the J term 
within the lowest order since J oc t 2 . This stimulates us to investigate in the following also separately the role 
of J term in Eq. (9), both through the MF term, Eq. (5), and the coupling Wk q , Eq. (3). 

NS results for Z^, ek are used as an input for the solution of the gap equation, Eq. (9), as presented in Fig. 2. 
For the same t'/t = -0.3 we calculate T c /t for c h = 0.12,0.17,0.22. Besides the result a) of Eq. (9) (full line 
in Fig. 2) we present also two alternatives: b) the solution of Eq. (9) without the MF term, and c) the result 
with mkq without the J term and omitted MF term. In the latter case, we used as input NS QP parameters, 
recalculated with correspondingly modified mkq- 

From Fig. 2 it is evident that the spin-fluctuation contribution is dominant over the MF term. When 
discussing the role of the J term in the coupling, Eq. (3), we note that in the most relevant region, i.e., along the 
AFM zone boundary mkQ = 2 J — 4£'cos 2 k x . Thus, for hole doped cuprates, t' < and J terms enhance each 
other in the coupling, and neglecting J in mk q reduces T c , although at the same time relevant is enhanced. 

The gap equation Eq. (9) leads to Ak with expected d x 2_ y 2 symmetry form Ak ~ Ao(cosfcj; — cosfc y )/2, 
with Ao(T = 0) ~ rjT c and r\ ~ 2.5. However, we observe that in the UD regime the effective coherence length 
£ <~ vf/Aq(T = 0) becomes very short. I.e., with vf taken as the average velocity over the region k at the 
antinodal part of the FS we get £ ranging from £ = 4.4 in the OD case, to £ = 1.3 in the UD example. In the 
latter case, SC pairs are quite local and the BCS-like approximation without phase fluctuations, Eqs. (9), (12), 
overestimates T c . Starting from this side, a more local approach would be desirable. Full numerical solution of 
Eq. (9) is due to strong momentum dependent coupling also momentum dependent beyond the simple <i-wave 
form. In Fig. 3(a) Ak s (T = 0) is presented throughout the Brillouin zone for underdoped Ch = 0.12 and coupling 
with emitted MF term as in (b) from Fig. 2 and for t'/t = —0.3. In the overdoped regime, Ch = 0.22, the 
momentum dependence of Ak s is less pronounced, Fig. 3(b), in accord with moderate momentum dependence 
of Zk, Fig. 1(b). 

As discussed above, we take here the dynamical spin susceptibility independent of temperature, the approx- 
imation justified for low temperature T ~ T c and solutions of the gap equation in the T < T c regime should be 
considered only qualitatively. In Fig. 4 full temperature dependence of Ak s (T) is presented for k = (ir, 0) for 
various c/j and t'/t = —0.3. An interesting observation is increasing ratio A^^ a (T — 0)/T c from OD to UD 
regime reflecting different relative relevance of J- and t'-terms in the coupling. 




Figure 3. (a) Zero temperature A ks obtained from Eq. (9) using coupling without MF term [as in Fig. 2, line (b)], for 
c h = 0.12 and t' ft = -0.3. (b) As in (a), but for c h = 0.22. 




Figure 4. Temperature dependence of Ak s (T) for k = (tt,0) with parameters and coupling corresponding to Fig. 3. 

Spin-fcrmion coupling mkQ is strongly i'-dependent and as expected also T c exhibits a pronounced dependence 
on t' It. In Fig. 5 we present results, as obtained for fixed OP Ch = 0.17, but different t'/t < 0, as relevant for 
hole-doped cuprates. 20 

4. DISCUSSION 

Let us first comment on the relevance of the present method and results to cuprates. Our starting point is the 
model, Eq. (1), where strong correlations are explicitly taken into account via the projected fermionic operators. 
In this respect the derivation crucially differs from the analysis of the phenomcnological spin-fermion model. 6 
Nevertheless, in the latter approach the resulting gap equation, Eq. (9), looks similar but involves a constant 
effective coupling. In contrast, our TOk q , Eq. (3), is evidently k, q- dependent. In particular, in the most relevant 
region, i.e., along the AFM zone boundary, mkQ depends only on t' and J, but not on t. This explains our 
central result novel within the spin- fluctuation scenario, i.e., a pronounced dependence of T c on t' which emerges 
directly via t' in the effective interaction in Eq. (12), and is consistent with the evidence from different families 
of cuprates. 20 Similar trend is obtained within the same model by the variational approach. 11 One can give 
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Figure 5. T c /t vs. —t'/t for fixed 'optimum' doping c h = 0.17 and different versions of Eq. (9), as in Fig. 2. 

a plausible explanation of this effect. In contrast to NN hopping t, the NNN t' represents the hopping within 
the same AFM sublattice, consequently in a double unit cell fermions couple directly to low-frequency AFM 
paramagnons, analogous to the case of FM fluctuations generating superfluidity in 3 He. 22 

It is instructive to find an approximate BCS-likc formula which simulates our results. The latter involves 
the characteristic cut-off energy Tq, while other relevant quantities are the electron density of states A/o and 
Z m being the minimum on the FS (at the antinodal point). Then, we get a reasonable fit to our numerical 
results with the expression, 

T c « ir Q e-W oV "'\ (12) 

where the effective interaction is given by V e // = 3Z m (2J — 4i') 2 XQ- Our numerical analysis suggests that the 
main i'-dcpcndcncc of T c originates in the coupling rn-kqi n °t in J\foZ m , while the main c^-dependence comes 
from xq and Tq. Then, Eq. (12) implies that optimum doping, where T c reaches maximum, increases with 
—t'/t. For parameters used in Fig. 1, e.g., c opt = 0.13 + 0.12(— t'/t). 

In this analysis we do not extend our input data outside the doping range 0.12 < Ch < 0.22. Nevertheless, 
we can discuss on the basis of Eq. (12) the variation T c (ch) elsewhere. Towards the undoped AFM also the spin 
fluctuation scale should vanish Tq — ► and consequently T c (ch — > 0) — > 0. On the OD side, xq and therefore 
V e ff should decrease with doping, leading again to fast reduction of T c (ch). 

It is evident from our analysis, that actual values of T c are quite sensitive to input parameters and NS prop- 
erties. Since we employ the lowest-order mode-coupling approximation in a regime without a small parameter, 
one can expect only a qualitatively correct behavior. Still, calculated T c are in a reasonable range of values in 
cuprates. We also note that rather modest 'optimum' T c value within presented spin-fluctuation scenario emerge 
due to two competing effects in Eqs. (9), (12): large rhkq and xq enhance pairing, while at the same time through 
a reduced Z^ and cutoff Tq they limit T c . 

It should also be noted that in the UD regime we are dealing with the strong coupling SC. Namely, we observe 
that AfoV e ff shows a pronounced increase at low doping mainly due to large xq- Then it follows from Eq. (12) 
that T c is limited and determined by Tq. At the same time, INS experiments 19 reveal that in the UD cuprates 
the resonant peak at w ~ u r takes the dominant part of intensity of q ~ Q mode which becomes underdamped 
possibly even for T > T c . Thus it is tempting to relate Tq to cu r (for more extensive discussion see Ref. 21) and 
in the UD regime to claim T c <~ Civ r , indeed observed in cuprates 19 with C <~ 0.26. However, additional work is 
needed to accommodate properly an underdamped mode in our analysis. 
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Figure 6. T c vs. ch with aF — (ch — co)i, b — 2at and c = 0. Dashed line connects the positions of Tmax = (c op t — co)i/e. 

Finally, let us conclude with a simplified qualitative picture emerging from our numerical analysis. As 
discussed above, full numerical treatment reproduces doping- and t'-dependence consistent with cuprates. NNN 
hoping matrix element t' enters the formalism in two ways. Firstly, through the dispersion of quasiparticlcs and 
corresponding renormalisation in NS leading to moderate increase of the density of states. Contrary to some 
other previous works, this effect is, surprisingly, not strongly doping dependent. More important is the effect 
of the spin-fermion coupling, where t' couples to spin susceptibility. The main ^'-dependence of T c comes from 
this part. The doping dependence emerges mainly due to strong spin susceptibility variation with c/j. Taking 
temperature independent Tq and xq oc Tq 1 together with the assumption of independence of the density of 
states one/,, a simple qualitative formula for T c emerges, 



As an example we take aT = (c/, — c )i, b = 2at and we neglect spin exchange in the coupling, c = 0. In 
Fig. 5 doping dependence of such T c is presented for various t' /t. In spite of severe approximations used in the 
derivation of Eq. (13), this result reflects the main features found also in the full numerical treatment: i) T c is 
maximum in OP and decreasing in UD,OD regime, ii) larger maximum T c at larger —t', and iii) optimum c opt 
increases with — t'. In the UD regime T is small compared to the relevant temperature scale and therefore for 
Ch ~ Co the use of a more precise form of T(ch,T) would be necessary. Maximum possible T c based on Eq. (13) 
is given with T max = 0.37ar, where bT = (c — t') 2 and therefore c opt = c + a(c — t') 2 / (bt). 

Such estimates are consistent with our numerical results and could be relevant also for real cuprate materials. 
Maximum possible transition temperature is then given with approximate formula 



T c = ar e - br /( c -*') 2 . 



(13) 



T max w 2- M Z m {2J - U') 2 XQ T Q . 



(14) 
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